Heteronuclear fermionic superfluids with spin degrees of freedom 
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We present a theory of spinor superfluidity in a two-species heteronuclear ultracold fermionic 
atomic gas consisting of arbitrary half-integer spin and one-half spin atoms. In particular, we focus 
on the magnetism of the superfluid phase and determine the possible phases in the absence of 
a magnetic field. Our work demonstrates similarities between heteronuclear fermionic superfluids 
and spinor Bose-Einstein condensates at the mean-fleld level. Possible experimental situations are 
discussed. 
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I. INTRODUCTION. 

Fermionic superfluidity is currently one of the most ac- 
tive research topics in the field of ultracold atomic gases. 
Fermionic superfluids have been realized and the BEC- 
BCS crossover regime has been explored for both bal- 
anced and imbalanced systems [H, 0, d, 13, H S S H, @ ■ 
Another line of research that has attracted a growing 
interest concerns mixtures of different atomic species 
[13, mi, [13, El, such as two-species Fermi-Fermi mix- 
tures. Fermionic heteronuclear gases have currently been 
inve stig ated both experimentally [T^j and theoretically 
[l^ llTl [isl . [TqI . Recently, the first quantum degenerate 
two-species Fermi-Fermi mixture has been realized [20| . 

In this paper we theoretically investigate superfiuid 
properties of a heteronuclear Fermi-Fermi mixture with 
spin degrees of freedom. Two-species Fermi mixtures dif- 
fer from the single-species gas in such a fundamental way 
that at low temperatures atoms in the same spin state can 
interact with each other through s-wave scattering. For 
identical fermions, however, this is prohibited due to the 
Pauli principle and consequently superfiuidity of a single 
species of atoms occurs in the spin-singlet state for the s- 
wave channel. Specifically, for the case where both atoms 
have spin one-half this means that both spin-singlet and 
spin-triplet states are allowed. Theoretical investigations 
of superfluidity of heteronuclear Fermi gases have so far 
not explored the physics of the (hyperfine) spin degrees of 
freedom of the atoms [l^] ■ While pairing in single-species 
fermionic superfluids with arbitrary spin has been stud- 
ied [ll], [l^] , no attempt to investigate the combined ef- 
fects of those two, i.e., heteronuclear spinor superfluidity, 
has been made to the best of our knowledge. This is the 
subject we address in this work. In particular, our work 
demonstrates a strong similarity between heteronuclear 
superfluids and spinor BECs at the mean-field level and 
suggests the existence of novel many-body states. 

Experimentally, a particularly interesting candidate 
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to realize a heteronuclear spinor superfluid is an iso- 
tope mixture of the rare-earth element ytterbium (Yb), 
which has two stable fermionic isotopes with nuclear spin 
7=1/2 (I'^^Yb) and 7 = 5/2 ("^yb), respectively, and 
electronic spin 5 = 0. Both isotopes have already been 
trapped optically and recently the ^^'^Yb gas was cooled 
to quantum degeneracy 123]. 

This paper is organized as follows. We formulate our 
problem in Sec. Ill Al and develop the Bogoliubov and 
Ginzburg-Landau theories for our system in Sees. Ill Bl 
and III CI respectively. We apply the theory for specific 
cases in Sec. IIIII and discuss possible experimental re- 
alization of the gases considered here in Sec. IIVI . We 
summarize and conclude this paper in Sec. |Vl Some 
algebraic manipulations to derive the Bogoliubov exci- 
tation spectrum of a heteronuclear superfluid have been 
relegated to the Appendix to avoid disgressing from the 
main subject. 



II. THEORY OF HETERONUCLEAR 
FERMIONIC SUPERFLUIDITY 



A. Formulation of the problem 

We consider an optically trapped Fermi-Fermi mixture 
with spins > 1/2 and = 1/2, where we distinguish 
the two atomic species with suffixes 4> and x- We assume 
that the temperature is sufficiently low so that we have 
only to consider s-wave interactions between the atoms. 
In general, there are interspecies and intraspecies inter- 
actions, and experimentally it is possible to tune the in- 
terspecies interactions independently of the intraspecies 
interactions. In this paper we will ignore the intraspecies 
interactions to focus on the bare essentials of this sys- 
tem. The total spin / of the interacting atoms with spin 
> 1/2 and = 1/2 is given by = U± 1/2. In the 
absence of an external magnetic field, the total Hamilto- 
nian H — Hib -I- is the sum of one-body part 77i_b and 
interaction V. The general form of the interaction can 
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be represented as 



with P 



V^J2 Vf{x-x')Vf, 



(1) 



/0 ± 1/2 that 
1 



Vf± = 
^ 2/H 



(5) 



where the projection operators V f project pairs of atoms 
onto pairs of the total spin / channel and Vf{x—x') is the 
interaction potential between a pair of atoms with total 
spin /. For s-wave interactions, the interaction potential 
can be approximated by a pseudopotential Vf{x — x') = 
gfS{x - x'). 

In second quantization the one-body Hamiltonian is 
given by 



The interaction Hamiltonian can then be constructed as 



V = [gf+Vf++gf-Vf-]d{x-x' 



+ 
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5{x - x') • /0. 



(6) 



dx 



The corresponding second-quantized expression is 



dx E ^i(^) 

CT=±l/2 



' 2Mv 



dx 



where (i/io-) and xt {xa) are the creation (annihilation) 
operators of atoms of spin and spin 1/2, respectively, 
in the magnetic sublevel a. 

The terms V^^^{x) and V^^^{x) describe state- 
dependent external potentials for the (j) and x atoms, 
respectively, fi^ and fi^ are the chemical potentials, and 
and are the masses of the and x atoms, re- 
spectively. The projection operators can be expressed in 
second quantisation as 



(3) 



where 

^f,m (x) 



(7, a' 



(/,n^l/0,cr; \/2,(T')(l)a{x)Xa'{x) 



is the annihilation operator of a pair of </> and x atoms 
with total spin / and total magnetic quantum number 
m, with {f,m\f^,a; 1/2, a') being a Clebsch-Gordan co- 
efficient. The interaction Hamiltonian is then given by 



V = 



2 

f.fn 



dx gfA)^^{x)Aj^^{x). 



(4) 



Since we are interested in the superfluid phase we assume 
the interaction between the atoms to be attractive, i.e., 
5/ < 0.^ ^ 

The interaction Hamiltonian in Eq. ([T]) can be rewrit- 
ten in a physically suggestive manner. It follows from 
the completeness relation 1 = 'Pf + + Vf- and from the 
total spin squared 



(/x + Uf = 



f=f+J- 



{9f+ -9f-) 



.9/ 



^F{x) ■ *'F(a;) 
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2/H 



9f 



2/^ 



(2) 



h^{x)n^{x) :, 
(7) 



where :: denotes normal ordering, h^{x) and fiy^ix) are 
the total densities of the (j) and x atoms given by 



and 



(8) 



(9) 



and ^F{x) and '^F{x) denote the spin density vectors 
whose components are given by 



and 



^F'{x) 



a, a' 



(10) 



(11) 



respectively, where \^F'^\„„t and [•*^-F'*]o-cr' are the matrix 
elements of the i — x,y, z components of spin matrix 
vectors "^F and ^F. The sign of {gj:+ — gf-) specifies 
which of the two spin states /+ and is energetically 
favorable and thus determines if the total angular mo- 
mentum of the Cooper pairs is /"*" or /~. The generic 
phase diagram that results is sketched in Fig. [TJ At high 
temperatures, the system is in the normal state and the 
total angular momentum is not fixed. When we lower the 
temperature for a fixed value gf- and \gf+\ > 1.9/- |, 
we expect a phase transition to a BCS superfluid whose 
total angular momentum of each Cooper pair is and 
the transition temperature is determined hy gf+. On the 
other hand, when \gf+ \ < \9f-\, we expect a BCS super- 
fluid with total angular momentum /~. By changing the 
ratio gf+/9f- we can switch between the two types of 
superfluidity. 
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with (F)^ ^, being the spin-/ matrix, and A/^m is given 
inEq. (fljj). 

In the absence of the spin degrees of freedom the Bo- 
goUubov modes in Eqs. ((M]) and describe quasi- 
particle excitations of fermions with mass and population 
imbalance [H, [13, 13 • To determine the spin structure 
we next derive and solve the gap equation for our system. 
After the Bogoliubov transformation the Hamiltonian is 
diagonalized as 



FIG. 1: Generic phase diagram as a function of temperature 
T and coupling constant gj:+ . ( /+) denotes the phase 
in which the spin angular momentum of a constituent Cooper 
pair is /" ( /+). 



B. Mean- field theory of heteronuclear superfluidity 

The phases we found in the previous subsection possess 
the internal spin degrees of freedom. To investigate the 
magnetic properties of the superfluid phases we make use 
of mean-field theory. The order parameter is defined by 



cr,k 



<T,k 



A/,™ (a;) = 5/(1/,™ (a;)). 



(12) 



Applying the Gorkov decoupling, we decompose the in- 
teraction Hamiltonian as 



|A/,™P 
2g/ ' 

(18) 

where r)cr,k and Co-.k are the Bogoliubov-mode operators 
that reduce continuously to (pa-^k and Xo-,k, respectively, as 
I A| 0. There are additional contributions to the above 
Hamiltonian arising from the degenerate unpaired modes 
with ideal gas dispersion ^^.k = Ti'k'^ /2M^ — /i^; these 
terms, however, do not depend on the energy gap and 
are therefore not relevant for calculating the gap equa- 
tion. The gap equation follows from the thermodynamic 
potential 17, which is determined from the partition func- 
tion Z = e~^^ = Tre~'^^. Using the Hamiltonian in Eq. 
(fTS)) . we obtain 



9f 



_ _ A/,„,p 

x,k '"^xM) 2gf 



(13) /3 
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In a homogeneous system, the resulting mean-field 
Hamiltonian can be diagonalized by means of a Bogoli- 
ubov transformation. We find that the Bogoliubov dis- 
persion relations are given by (see Appendix[X]for deriva- 
tions) 



(19) 



The gap equations that follow from dft/dA*^^^ = are 
given by 



A 



and 



(14) 
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T: 
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f l=±l k 

iV(nt^' )+7V(?ic.t,k)-l 



2 e0,k + ^xM 



2 (?it^0_k + 



xM> "Ma . ' ^l^'^l 



X.k 



2/+ dA 



(20) 



where = e^^k - ^^a, e^.k = ?i^A;^/2M„ (a = 0, x), 
and 



^^±-^y^|Ap±^ + (e0,k + ex,k)'. 

(16) 

Here |Ap = Z]m|A/,mP, the spin vector Fa of the 
Cooper pairs is given by 

Fa = At . F . A = X A},„ (F),„_„, A/,„,, , (17) 



Here by way of renormalization we have replaced the bare 
coupling constants f// with the two-body T-matrix T/^, 
where T'^^ is related to 5/ by 



rp2B 



9f V 



(21) 



From an experimental point of view, the two-body T- 
matrix is directly related to the s-wave scattering length 
a/ with the total angular momentum / of the two 
colliding atoms via T^^ = 2'Kafh'^ /M^, where = 
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M^M^/{M^ + M^) is the reduced mass. If A/,„ 7^ 
we can divide both sides of the gap equation by Ay ™ 
and, as a result, the left-hand-side of the gap equa- 
tion (PO)) becomes independent of m. To be consistent, 
the right-hand-side should also be independent of to, so 
'^E^ must be independent of m. This imposes a 
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self-consistency relation for the order parameters and de- 
termines the possible superfluid phases. Before we inves- 
tigate concrete examples, we will first investigate the con- 
nection of heteronuclear fermionic superfluids with spinor 
BECs in the next section. 



C. Ginzburg-Landau theory 



spin can be maped onto a spinor BEC. For this pur- 
pose we formulate the Ginzburg-Landau theory for our 
system. We make use of the functional-integral formal- 
ism and as our starting point we take the action for the 
Hamiltonian Hib + V discussed in Sec. Ill Al The parti- 
tion function Z is expressed in terms of functional inte- 
grals over the fermionic fields (pa and Xc ^ 



Z 



d[4>]d[c^*]d[x]d[x' 



-S/h 



(22) 



In this subsection we show that a heteronuclear 
fermionic superfluid with Cooper pairs having a nonzero 

I 



where the action is given by 



sWA,x\x\ 



J dr J dx'^Xa(x,T) (^i^^ 



d 



2M^ ' X 



+ V^'^{x)-^i^]xA^,r) 



1 f^'^ f 



(23) 



The first two lines on the right-hand side describe free 
propagations of the fermion fields and the last line de- 
scribes the two-body interaction which involves pairing 
amplitudes 



Af^„i{x,T) = 

^{f^ Mf<l>, ^,Cr')<l>a{x, T)Xcr' {X, t). 
ct.ct' 

(24) 



We consider the homogeneous case here by setting the 
external potentials V°^g.{x) and V^^^{x) to be zero. 
We decouple the interaction by means of a Hubbard- 
Stratonovich transformation, i.e., we introduce a com- 
plex auxiliary field A/^„i(a;,r) that couples to the prod- 
uct of fields X;^,<T'(/''^I/0'O-; i,cr')(?!)^(a;,T)xa'(a;,r). 
The fourth-order term generated by the Hubbard- 
Stratonovich transformation exactly cancels the inter- 
action term in the action and the resulting action 
S[A* , A, (/)* , (j),x* ,x] depends only quadratically on the 



fermion fields <paix, r) and Xa-{x, t): 

S[A*,A,r,cb,x*,x] = - f^dr ^dxY,- 
•'" J m 



9f 



(25) 



where ^* stands for a set of fields 

** {x, t) = (^(f)}^ {x,t),... 01 j:^ (a;, t),Xi{x, t), X^Jx, t)^ 

and properties of the system are encapsulated in the 
Green's function matrix G which is expressed in terms of 
the noninteracting Green's function Go and self-energy 
S as 



(26) 



G-\x, t; x', t') = Go\x, t; a;', t') - S. 



(27) 



The noninteracting Green's function Go is the diagonal 
(2/0 + 3)x (2/0 + 3) matrix 



G-\x,t;x',t') = 
G^lix,T;x',T')l 



-G-Ux\t';x,t)1 



(28) 
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where -G^'^q{x' ,t';x,t)1 {a = (j),x) is a (2/^ + 1) x 
(2/0 + 1) matrix that satisfies 

= _i InS- - %^ - Ma) - x')5{t - t') . 



(29) 



The self-energy S is a (2/^ + 3) x (2/^ + 3) matrix 

(30) 

where V is the (2/^ + 1) x 2 matrix given by Va-.a-' = 
cr; 1/2, o-')A/^„i. Since the action (gS]) de- 
pends only quadratically on the fermion fields, we can 
integrate them out and obtain an effective action 



S"''^[A*,A] 



h/3 



dr / dx 



|A/,„(a;,r) 



9f 



The last term in Eq. pip can be expanded in powers of 
A by using 



G-^=Go-^-S = Go^(l-GoS), 



-;iTr[log(-G-i)] . 



(31) 



where the self-energy KTi is given by Eq. ((30|) , and there- 
fore 



;iTr[log(-G-i)] = -?iTr[log(-Goi)] 

+h V -Tr[(GoS)™] . (32) 
^ — ' m 

m— 1 



The second-order m = 2 term in the last term on the 
right-hand-side of Eq. ([5^ is given by 



-Tr 



(GoS)^ 



dr dx dr' / dx' / dr" / dx" / dr'" / dx'" 



tr [Go(a;, t; a;', r')S(a;', r'; a;", t")Go{x" , t"-x'", t"')T.{x'" , t'"; a;, t)] 



(33) 



where the trace operation tr[. . .] on the right-hand side 
means that we only take the sum of the diagonal ele- 
ments of the (2/0 3) x (2/0 + 3) matrix GqSGoS. In 
the following we assume A/^™ to be independent of space 
and (imaginary) time. This allows us to separate the 
imaginary-time and space integrations from the matrix 
trace. For equal chemical potentials and equal masses 
we have G~^\{x, r; x' , t') — G^^ix, r; x' , t'), and the in- 
tegrals over the Green's functions reduce to those of the 
standard BCS theory. 



The new contributions come from the matrix trace. 
For the second-order term we find after some straightfor- 
ward algebraic manipulations that 



To evaluate the fourth-order term, we make use of 
tr[yV\/V] = tr[V^FVFt] 

, Fa -Fa 



(35) 



Combining these with conventional BCS theory we find 
that to the fourth order the thermodynamic potential, 
for equal masses and chemical potentials, is given by 



n = a/(T)E|A;,„,|2 + /3 5]|A/,, 



\Fa\ 



tr[Wy] =tr[T/W] =E|A/,™|^ 



(34) 



(36) 



where af{T) = Ar(0) logr//T, /? 
7C(3)iV(0)/16(7rfcBT/)2, N{{)) = 2MkY/{2T:hf is 
the density of states at the Fermi energy, which is the 
same for both species since M = Ms = My [ij. The 
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critical temperature is given by 
Sep 



g7-2g-7r/2|fcFa/ 



(37) 



The form of the thermodynamic potential in Eq. (j36[) is 
essentially the same as that of the homogeneous spinor 
BECs at zero magnetic field. Here, however, since f3 > 
the ground state is unmagnetized, i.e.. Fa = 0. Physi- 
cally this result arises from the Pauli principle which for- 
bids the fermionic constituents of the composite bosons 
to occupy the same state. 

Moreover, we only find two kinds of fourth-order terms, 
i.e., the particle density and spin density, regardless of 
the value of the spin /. In contrast, for the case of spinor 
BECs there are / +1 different fourth-order terms. There- 
fore, all unmagnetized phases predicted in spin-/ BECs 
are degenerate in the present case. 



III. RESULTS 

As shown above, the heteronuclear fermionic supcrfluid 
has a Ginzburg-Landau free energy which is similar in 
form to that of a spinor BEC. We may thus expect sim- 
ilarities between the two systems insofar as mean-field 
theory is valid. To investigate this problem, we solve 
the gap equation (|20p for the case of two spin- 1/2 species 
whose total spin is either zero or one with the correspond- 
ing s-wave scattering length being oq or oi. The phase 
which the gas condenses into is determined by the s-wave 
scattering length with the largest amplitude. For f = 
we immediately conclude that Fa = 0. The Bogoliubov 
spectrum is the same as that of the conventional BCS 
theory. For total angular momentum / = 1 we have 



E 



A 



l.m I 



3|ei| 



(38) 



where the scalar quantity 9i = (— q + 
2Ai_iAi__i)/-\/3 is the same in form as the spin- 
singlet state pair amplitude of a spin-1 BEC [s^l- If all 
Ai^m are nonzero, we find from the gap equation that 
either 8i = or 

1 961 1 dOi 



Ai,„. aA* „ Ai,„, aA* , 



(39) 



holds for m, m' = 0, ±1. It follows from the latter condi- 
tion that |Ai_i| — |Ai__i|. Moreover, if we represent the 
order parameters as Ai „j = |Ai „i|e*^^ '", we find from 
the gap equation that 26i,o = + 6*1,-1 +7r. We also de- 
fine the angle /3 from the relation Ai_o = (|A|/-\/2) cos/?, 
where |Ap = X]ml^/,mP- Combining these, we find 
that 



(Ai,i,Ai,o,Ai,_i) 



lAL 



1 



^e*" sin/3,cos/3, j=e '"sin/3 



V2 



V2 



where a = Oi^i — 9i^q. This order parameter also de- 
scribes the polar phase of the spin-1 Bose gas. From 
the condition Oi = 0, we obtain the phase relation 
26*1,0 = 6*1,1 +6*1,-1 and the amplitude relation |Ai,oP = 
2|Ai,i||Ai,_i|. From these results we find that the ferro- 
magnetic order parameter is given by 



(Ai,i,Ai,o,Ai,_i) = 



e cos — , V 2 sm — cos — , e 
2 2 2 



^ . 2/3 
sm — 
2 



(41) 



(40) 



It can be shown that Eq. (|4T|) is indeed a solution to 
the gap equation (PO)) . Comparing the free energies of 
the polar and ferromagnetic phases we conclude that the 
ground state is polar. Because the order parameter of the 
polar phase has a Z2 symmetry, i.e., is invariant under 
the transformations a a + n and /3 — > tt — /3, the mean- 
field theory predicts a quantum phase transition when we 
change the scattering lengths from ai > ao to oi < ao, 
and vice versa. 

To find solutions of the gap equation for higher val- 
ues of the spin /0 > 1/2 it is convenient to make use 
of the one-to-one correspondence between the fermionic 
heteronuclear gas and the spinor BEC derived in the pre- 
vious section. The solutions for total spin 2 and spin 3 are 
of special interest because the correspond to the exper- 
imentally realized case of a mixture of ^^^Yb (spin 1/2) 
and "^Yb (spin 5/2). The spinor BECs with spin 2 and 3 
have been described in the literature 0, HI, [H, US] ■ and 
it is straightforward to verify that the possible ground- 
state solutions of the spin 2 and spin 3 BECs also satisfy 
the gap equation (|20|) . For example, for the case of the 
spin 2 BEC cyclic and polar phases are possible [l^, . 
These phases can be distinguished from each other by 
the value of 62 = {Al o - 2A2,i Aa^-i + 2A2,-2A2,2)/V5, 
which describes the formation of singlet pairs of spin-2 
atoms. For spin 3 the solutions are more complex and 
have been discussed by Santos and Pfau ^26.] and Diener 
and Ho [l^l . Also a complete classification of states and 
vortex excitations have been discussed in Refs. |28|. 12911. 



IV. POSSIBLE EXPERIMENTAL 
REALIZATIONS 

The critical temperature needed to reach the weak- 
coupling limit is beyond experimental reach and in order 
to increase the critical temperature to experimentally ac- 
cessible values we need to make use of optical Feshbach 
resonances [30] so that the effective interactions between 
the atoms are enhanced. Currently, optical Feshbach res- 
onances for alkali atoms [32,, i33j are not as effective as 
magnetic Feshbach resonances due to large atomic losses. 
It has been predicted [U, [11], however, that for ground- 
state alkaline-earth-metal atoms optical Feshbach reso- 
nances can be used to tune the scattering length over 
a wide range of values without suffering from the large 
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atomic losses. As mentioned in the introduction, the 
fermionic isotopes of Yb have already been trapped opti- 
cally [2^ , which makes them a prime candidate for realiz- 
ing a heteronuclear fermionic superfluid with spin degrees 
of freedom. At present, experimental efforts are under- 
way to find appropriate bound states in the long-lived 
(> lOsec) excited state potentials ^ So+^Pi=o 2 of Yb that 
can be used for optical Feshbach resonances [s^ . If these 
efforts are successful we expect that the ground-states of 
the ^^^Yb and ^^'^Yb mixture are the same as the unmag- 
netized phases of the spin 2 and 3 spinor BECs. More- 
over, in the absence of symmetry breaking perturbations, 
the excitation spectrum for the unmagnetized ground- 
states is given by Eqs. (fT4|) and (fTS]) with |Fa| = and 
has the same form as the standard BCS theory. For a 
realistic gas, however, symmetry-breaking terms such as 
magnetic fields and dipole-dipole or intraspecies interac- 
tions are expected to lift the degeneracy and we plan to 
investigate this in future research. 



SUMMARY AND OUTLOOK 



In summary, we have presented a theory of superflu- 
idity in heteronuclear Fermi-Fermi mixtures. We have 
shown that the heteronuclear fermionic superfluid has a 
Ginzburg-Landau free energy which is similar in form to 
that of a spinor BEG and, therefore, analogous mean-field 
ground states are expected for these two systems. This 
raises an important question concerning the nature of the 
many-body ground state in our case. For the bosonic case 
it has been shown that the exact many-body state is, in 
fact, not the polar state but a condensate of spin singlet 
pairs [H, m, \^ . Explicitly, in the absence of a magnetic 
field the exact bosonic many-body state can be written as 

((aj)2 - 2alaLi)/V3 |0), where aj„ creates a 

boson in hyperfine state m = 0, ±1 with zero linear mo- 
mentum and N is the number of bosons. To perform a 
similar many-body analysis for the heteronuclear Fermi 
superfluids is beyond the scope of the present paper and 
will be left for future research. It is clear, however, that 
if the similarity at the mean-field level extends to the 
many-body case, the exact ground state is expected to be 
a condensate of spin-singlet states involving a quartet of 
fermions. Quartet superfluidity is of interest in strongly 
interacting quantum liquids [4ll . |42| | and in the formation 
of two-pion states fi^ . In these works, however, only sin- 
gle species gases were considered for which a fourfold de- 
generacy and moderate coupling are indispensible. Our 
work demonstrates the possibillity of creating a quar- 
tet condensate using only spin- 1/2 particles because the 
Pauli principle does not prevent different species from 
occupying the same spin state. This unique feature of 
heteronuclear fermionic superfluidity merits further the- 
oretical and experimental investigations. 
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APPENDIX A: THE BOGOLIUBOV SPECTRUM 
OF A HETERONUCLEAR FERMIONIC 
SUPERFLUID 

Substituting the mean-field expression in Eq. (|13p into 
Eq. ([2]), we obtain the following Hamiltonian for a ho- 
mogeneous system: 



/ 



H 



E 

k 



14, 



Xl/2 
V-I/2J 



(e0,k-A)i V 



A)l 



9f 



xl/2 



(Al) 



where we define the (2/0 + 1) x 2 matrix 



V 



( y,M2 



^J-l.-l/2 



v^-,.1/2 vi,-_i/2y 



(A2) 



with matrix elements V'o-.o-' — 

I]m(/''^b'<^; The eigenvalue prob- 

lem associated with this Hamiltonian can thus be 
written in a compact form as 



(^k - A')l V 
1/t (-Ck-A')! 



0. 



(A3) 



Here, we have introduced the variable ^k = (^,^,k+Cx,k)/2 
and the shifted eigenvalue A' = A — (^^^k — ^x-k)/2- The 
vector u has 2j + 1 components and v has two compo- 
nents. The highest order of A' in the eigenvalue eqation is 
^'2j+3 ^ determines the maximum number of modes 

in the system. After working out the multiplications, we 
obtain the following coupled equations 



(a-A')lw + l^-i' = 



(A4) 
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and 



- A')li' + V'' ■u = 



(A5) 



We solve Eq. (jA4p formally for u and substitute the 
solution into Eq. (|A5|) to obtain the following equation 
for v: 

(-a - A')l + ■ ^^^i ■V^v = 0. (A6) 

The eigenvalues of this equation can be obtained by solv- 
ing the characteristic equation, 



Det [(Ck - A')(-a - A')l + V^V] - 0. 
Similarly, we can obtain for u that 

(a - A')i + V ■ • i^U w = 0. 



(A7) 



(A8) 



The eigenvalues for this can be obtained by solving the 
characteristic equation 



Det [(a - A')(-a - A')l + VV^ = 0. 



(A9) 



While Eq. (KT} is a 2 x 2 eigenvalue equation, Eq. (|A9l) 
is a (2/0 + 1) X (2/0 + 1) eigenvalue equation which has 
2/0 — 1 degenerate solutions A = ±a in addition to the 



nontrivial solutions. The nontrivial solutions follow from 
Eq. dUl). We have, 

y^y ^ f Er ^;,i/2^jm/2 "(^;,i/2^/,-i/2 
\T,f -1/2^/, 1/2 "^^y, -1/2^/, -1/2 

(AlO) 

Using this, we obtain the eigenvalues A from Eq. (|A7p . 



CkA^ , 1 



±y(^Yj + ^tr[ytv-] ± i^(tr[yty])2 _ 4Det[l/tF]. 

(All) 

Next, we express the trace and determinant terms ap- 
pearing in the above eigenvalues in terms of physical 
quantities. Firstly, combining Eq. (jAlOp with Eq. (|A2p . 
we obtain 

tr[Fty] ^^|y^.,^,,|2^^|A^,„|2 (ai2) 



The other term is given by 

(tr [V^V]Y -4Dct[V^V]=^-^, 



(A13) 



where -Fa is given by Eq. pT)) . 
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